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^Nj Abstract 

^ We describe a system of stochastic differential equations (SDEs) -which model the interaction bet-ween 

O processive molecular motors, such as kinesin and dynein, and the biomolecular cargo they tow as part 

of microtubule-based intracellular transport. We show that the classical experimental environment fits 
within a parameter regime which is qualitatively distinct from conditions one expects to find in living 
^ ^ cells. Through an asymptotic analysis of our system of SDEs, we develop a means for applying in vitro 

observations of the nonlinear response by motors to forces induced on the attached cargo to make ana- 
q/^ lytical predictions for two parameter regimes that have thus far eluded direct experimental observation: 

1) highly viscous in vivo transport and 2) dynamics when multiple identical motors are attached to the 
• cargo and microtubule. 

g 1 Introduction 

Critical to the proper functioning of a biological cell is the efficient internal transport of organelles and other 
intracellular cargo that collectively form the basis of the cellular infrastructure [25,. In eukaryotic cells, 
-which are characteristic of all complex organisms, these biological materials are generally assembled near the 
nucleus and are then packaged into membrane-bound vesicles to be distributed through the cytoplasm to 
\^ appropriate locations throughout the cell. For larger vesicles and organelles (on the order of 1 fim) , diffusion 

is too slow for efhcient transport, and inadequate for proper spatial distribution of important organelles and 
' molecules throughout the cell. Cells have therefore evolved an intricate transport apparatus consisting of a 

7-H cytoskeletal net-work of thin filaments called microtubules and proteins called processive molecular motors 

T-H that move in a directed fashion along the microtubules -while generating sufficient force to to-w vesicles, -which 

-we -will generally call cargo, at a considerable rate (as high as 1 ^m/s |48l 1571 
^ The class of motor proteins we consider in this -work are kinesins, -which are typically responsible for 

transport to-ward the cell periphery. T-wo decades of biochemical and biophysical investigations have lead to 
a -working model for ho-w motor proteins move along microtubules f49| : Kinesins are dimeric molecules con- 
sisting of two heads that attach to a microtubule, a coiled-coil tether that joins the two chains, and a cargo 
binding tail (Figure [l| Each head contains both an ATP and a microtubule binding site, and these mo- 
tors "walk" along the microtubule track by chemical coordination of ATP hydrolysis cycles such that the cycles 
remain out of phase and at least one head remains bound to the microtubule. We say the motors are "proces- 
sive" because they take many steps during each encounter with a microtubule. All told, we have substantial 
information about the operation and properties of individual molecular motors bound to cargo through 
work done by biologists 13, ^ i2Sl [21] , physicists |571 [Ml 12] , and mathematicians Pl [Ml [55} [51 [5^ [^ . 
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However, translating these experimental observations and the as- 
sociated theory for single motor and cargo systems into a predictive 
model for in vivo intracellular transport has been constrained by 
several practical limitations. First, the two primary experimental 
techniques, described in Section [2] are conducted in fluid environ- 
ments that have substantially lower viscosity than what is expected 
in the cytoplasm [l^. Second, it is expected that the in vivo trans- 
port of a given cargo will involve multiple motors, possibly of dif- 
ferent types, that attach and detach from microtubules dynamically. 
Indeed, in their pioneering Markov chain model |45lll51H7] . Miiller, 
Klumpp and Lipowksy indicate the remarkably complex behavior 
that can result from a cargo responding to molecular motors ar- 
ranged in a "tug-of-war" configuration. However, with one notable 
exception where motors were artificially bound to the cargo in pairs 
[28| , it is experimentally difficult to infer precisely how many motors 
of each type are bound to both the cargo and a microtubule at the same time. In recent years, there- 
fore, theoretical work and numerical simulations of multiple-motor-cargo systems have both played a strong 
role in generating predictions that have been later verified, or modified as appropriate, by experimental 

findings [311 ig ini Els 153 [MSZl [HI IIHJ 137115115^ 
1.1 Outline of the Paper 

Our goal is to contribute to these efforts by constructing a stochastic dynamical system that is both analyt- 
ically tractable and rich enough to model the diverse conditions relevant for microtubule-based intracellular 
transport. Our system of stochastic differential equations (SDEs) very naturally incorporates the follow- 
ing qualitative features of motor dynamics: a linear or nonlinear spring model for the motor-cargo linkage, 
thermal fluctuations that influence the cargo, a nonlinear response by the motor to cargo forces, nuanced 
changes in behavior due to motor configuration, and fiuctuations in the motor position arising from the 
stochastic timing of the various mechanical and chemical events that occur during motor stepping. These 
features were absent in the original MKL models [451 146| . but are now common, not always in combination, 
in recent simulation-based studies ^ ,721 [311 ES] . We find that they substantially affect the predictions of 
the model. 

The framework we develop is quite general, but we restrict the model considerably in order to perform 
specific rigorous analysis. In particular, we study the transport of thermally fluctuating cargo by identical 
cooperative motors in the absence of motor binding and unbinding. This leads to certain results that are 
in contrast to prior studies, but when this happens, our analysis suggests that such properties may be a 
consequence of attachment/detachment dynamics which warrant separate study (see for example, |45|). Our 
purpose in neglecting binding and unbinding dynamics in the present study is not to dispute their relevance, 
but to provide a point of reference for understanding how the various properties of the motors and cargo 
affect their collective behavior. 

The core of our approach lies in Section [3] where we perform a dimensional analysis of the system of 
SDEs with physically relevant parameters. We find that there is a crossover in qualitative behavior between 
the diffusion-dominated, low viscosity in vitro environment in which many experiments are conducted, and 
the high viscosity environment in vivo. This result follows from the identification of critical nondimensional 
parameter groups summarized in Table [2j In the low viscosity regime, the rapidly fluctuating cargo achieves 
a quasi-stationary distribution with respect to the slower moving motors; in Section [3] we introduce the 
proper stochastic averaging techniques to calculate the influence of the cargo on the motors, see ([s]). In 
Section |4.1| we demonstrate that this separation of time scales results in a motor-cargo system velocity that 
is insensitive to order-of-magnitude changes in the viscosity of the fluid environment. Indeed, this has been 
observed in vitro by Gross et al [16] and in silico by Kunwar and Mogilner |36| . This framework also allows 
us to easily predict the result of perturbing various other parameters. For example, with the development 




Figure 1: A cartoon model of intra- 
cellular transport. Intracellular cargo 
may be attached to microtubules by 
one or more molecular motors which 
individually tend to walk toward ei- 
ther the plus-end or minus-end of a 
microtubule. 
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of synthetically modified motors |581 [751 15T] it may soon be possible to tune the spring constant of the tail 
of the motor which binds to the cargo. From equation ( 13 1 we observe that as the spring constant increases, 
the motor-cargo system velocity will decrease. 

Recently there has been some discussion about how to translate inferences made from low viscosity in 
vitro experiments into predictions more like the conditions of a biological cell [TFl IMl H5| . The cytoplasm 
is highly viscous, and likely to have some viscoelastic properties. It is beyond the scope of this paper to 
address diffusion in a viscoelastic environment, so we adopt the same approximation for cytoplasmic viscosity 
as presented by Mitchell & Lee |12 • In Section |4] we attempt to establish an equivalence between the drag 
force felt by the motor via the cargo at high viscosity with forces artificially introduced by experimenters 
in a low viscosity regime. This leads to a general formula (16 1 for calculating the velocity of a single motor 
with cargo as a function of solvent viscosity, and we provide a "back-of-the-envelope" formulation in equation 
(18 1. While these formulas provides a convenient interpolation between in vitro and in vivo velocities, our 
model does predict that the force required to stall a motor increases as the viscosity of the environment 
increases. This counterintuitive result is not in agreement with the experimental work of Shubeita et al |64| . 
Resolving this disagreement merits a more detailed analysis of the near-stalled motor mechanics, which we 
hope to address in future work. 

Having established these results for a single motor interacting with cargo, we proceed in Section |5.1| to 
analyze the dynamics when two identical motors are attached to the cargo. There are two primary mecha- 
nisms by which one may expect multiple motors to enhance transport. First, cargo bound to microtubules 
by multiple motors are able to stay attached, and hence stay in transport mode, longer. Second, because the 
cargo load is distributed among the motors, each individual motor carries less load and should be able to 
process along the microtubule more quickly. However, it has recently been observed that under light cargo 
loads, transport with multiple motors may be slower than with one |37L [7^ [5S] . 

We consider the latter phenomenon by way of two averaged models. In Section 5.1 we study an "instant 
relaxation" model where the cargo immediately moves to a position centered among the motors, a simplifying 
assumption also made in |72| . This is essentially equivalent to assuming the time scale of the cargo dynamics 
is fast compared to the motors (which we show in Section |3] is typically true in vitro, but not necessarily in 
vivo), and that thermal fiuctuations of the cargo can be neglected (this latter assumption does not have a 
similar quantitative justification). We compute transport properties of the motor-cargo system as a function 
of external force. In this simplified setting, we show that whether or not two motors are slower than one, 
in fact, depends on the shape of the force-velocity curve through its influence on the spatial configurational 
dynamics of the motors. This result is formalized in Theorem |5.1[ We also demonstrate that the force 
required to stall a two-motor system is more than twice that of a one-motor system. This superadditivity of 
stall force is in contrast to prior studies, suggesting that previously observed suhadditivity [36' may again be 
a consequence of attachment/detachment dynamics. 



2 The mathematical model 
2.1 Experimental Considerations 

There are two primary experimental techniques for studying how teams of motors respond to a given load. In 
microtubule gliding assays [14,, the tails of kinesin, which are typically bound to cargo, are instead adhered 
to a glass plate. A microtubule is placed on the motor-coated plate; by binding to and taking steps on 
the microtubule, the kinesin motors generate sufficient force to transport the microtubule, whose average 
velocity is observed. Optical tweezers [52 are used to apply an opposing force to the trailing end of the 
microtubule while its resulting velocity is monitored. In a second type of experiment, optical tweezers are 
used to track latex beads, serving as in vitro cargo, which bind to an unspecified number of motors and are 
transported along individual microtubules. Through an active feedback system, the tweezers are capable 
of applying approximately constant force while tracking the beads. We will focus primarily on this second 
experimental setup. 

In nearly all of the previous theoretical efforts, experimental data for single motor-cargo systems has 
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been used to build models for the behavior of motor-cargo systems with multiple motors. In our model as 
well, experimental data for single-motor cargo systems plays a vital role. We briefly summarize the most 
important qualitative properties of single motor and cargo systems: The motors take discrete steps along 
the microtubule - for kinesin-1, approximately 8 nm [7]. The stepping rate depends on the local ATP 
concentration [57] and on external forces applied to the cargo |671 15] . The coiled-coil tail between the motor 
heads and the cargo is a semi-flexible elastic structure |7]- Eventually, the motors will unbind from the 
microtubules, the cargo, or both. The stronger the opposing force applied to the cargo, the shorter the run 
length before such detachment occurs [591. 

In terms of the motion of motors along the microtubule, many existing models employ discrete state-space 
Markov chain methods [721 [SJ [SSJ [TSJ [S]. For example, in [72 , motors undergo random walks on a spatial 
lattice, with hopping rates biased according to the strain in their tails connecting to the cargo; the cargo, in 
turn, is assumed to instantaneously relax to a minimal energy position conditioned on the current position 
of the motors. A lattice is well-suited for the discrete-state stepping process of the motors and lends itself to 
incorporating both the biochemical transitions and the continuous spatial dynamics of the unbound head and 
cargo. However, the continuum model we choose is a natural coarse-graining obtained through systematic 
asymptotic reduction procedures |70l ITTl [T^ IITI [571 [?7l f5T| . In particular, starting from a more detailed 
dynamical description for the molecular motor, formulas for the drift and diffusion coefficients appearing in 
Eq. ([ij can be computed through Markov chain calculations |70l [7T1 [T^ , renewal theory [411 1571 l?fl [25] , or 
homogenization theory |51| . 

Rather than exploring the discrete-continuous connection in detail here, however, we build our model at 
the mesoscale, based on single-motor force-velocity and force-diffusivity relations. The qualitative features 
of these relations are informed by the robust literature of experimental measurements [671 1521 [521 [55] . 
Furthermore, while the models of [371 [Ml [251 [551 172] do have spatial detail and incorporate force-velocity 
relationships comparable to our model, our formulation permits a more flexible and faithful rendering of 
the random diffusive component of the motor dynamics. We find that this diffusive component, along 
with the shape of the force- velocity curve, plays a crucial role in Proposition |5.2[ where we show that the 
model exhibits a superadditive growth of the stall force for a collection of motors. Our choice of modeling 
framework renders asymptotic analysis possible; we identify which parameters are physically small and use 
stochastic averaging techniques to derive a simplified, lower-dimensional system of equations whose solution 
approximates the dynamics of the original model. In particular, for the case of TV = 2 motors, we use this 
systematic dimension reduction to obtain a completely analytical relationship between the dynamics of the 
motor-cargo complex and the model for a single motor. 

2.2 An SDE modeling framework 

The location of a molecular motor along a microtubule as a function of time t is described by a one- 
dimensional spatial coordinate Xi(t), indexed by the motor i e {1, . . . ,N}. In particular, we coarse-grain 
over details of the biochemical and conformational states of the motor, and we consider length scales (on 
the order of hundreds of nanometers) sufficiently large that the motor position can be reasonably described 
by a continuous spatial variable rather than in terms of discrete steps of a few nanometers. Because we 
consider the case of N cooperative and identical motors, each motor has identical dynamics. We represent 
the position of the cargo by a one-dimensional spatial coordinate Z{t) along the microtubule, neglecting 
transverse fluctuations. Korn el al |34| numerically explore a similar spatial continuum model that allows a 
three-dimensional representation for the cargo, but since our purpose here is to apply stochastic averaging 
techniques to derive analytical coarse-grained descriptions, we choose to illustrate this process in the simplest 
case. 

In this model, the dynamics of each motor and the cargo to which it is bound are described by the 
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following system of stochastic differential equations 

dX,{t) = vg{F{X,{t) - Z{t))/F,)dt + (Th{F{X,{t) - Z{t))/F.,) dW,{t) 



N 



ldZ{t)^ \j2F{X,{t),Z{t)) -e\At+ ^2kBT-fdW,{t), (1) 

i=l 

where is the stall force of a motor; v is the mean speed at which a motor moves along the track when 
no force is being applied; |cr^ is the effective diffusivity of a motor along the microtubule when no force is 
applied; 9 is the magnitude of any laser trap force applied to the cargo; T is the absolute temperature of the 
system; and fee is Boltzmann's constant. 

For the friction constant of the cargo 7, we use the Stokes-Einstein drag law for spherical particles, 
7 = Girar], where a is the effective radius of the cargo and rj is the dynamic viscosity of the fluid environment. 
There are in fact several appropriate choices for 7 when one considers the various types and shapes of cargo 
found in vivo, see Mitchell & Lee [15] for a thorough discussion. Our default simulations will be based on 
spherical 500 nm particles in a solvent with the viscosity of water at 300 K. This is appropriate for describing 
a typical in vitro experimental setup. The cytoplasmic environment, however, is thought to have a viscosity 
as much as 600 times that of water , in Sections |4] and [s] we will study our model across a broad range 
of values of viscosity. Rather than explicitly representing the results separately in terms of cargo size a and 
viscosity rj, which then would require a constant appeal to a spherical cargo assumption, we will present all 
results in terms of the more broadly applicable lumped friction constant 7. 

The nondimensional functions g and h are rescaled force-dependent drift and diffusion functions for the 
motors. The drift determines the instantaneous expected velocity of the motor and the diffusion function 
determines the local diffusivity of the process. Here we will simply take the diffusion as independent of the 
applied force {h = 1), and discuss the drift function in greater detail in Section [ZSj Typical values for the 
dimensional constants appearing in ([T]) are as listed in Table 2.2 for kinesin-1 



The functions {Wi{t)}fLi, Wz{t) are independent standard Brownian motions: that is, each of {Wi{t)}^i, Wz{t) 
is a continuous, Gaussian stochastic process with independent increments, zero mean, and variance equal to 
t; and Wz{t) are independent. There is a noteworthy distinction between modeling justification 

of the process Wz{t) and the system of processes {Wi{t)}^^. The SDE for Z is the overdamped Langevin 
equation for a diffusing spherical particle, subject to external forces that appear in the drift term. Were it 
not for these forces, Wz{t) would represent the Brownian motion of a free particle. By contrast, the processes 
{Wi{t)}^i represent the deviations from the mean of the motor dynamics that arise from fluctuations in 
chemical reaction wait times (and other discrete random events) associated with stepping. For this reason, 
the SDE for the cargo should satisfy the fluctuation-dissipation relationship via the drag parameter 7, while 
the motor equations need not. 

The function F{r) describes the force exerted by the cargo on a motor through its connecting tail, in the 
direction against its natural forward motion, where r is the signed separation between the motor and cargo 
positions. As in the literature, we use spring-like models for this force law. The simplest Hookean model 

F{r) ^ KT (2) 

with constant spring coefficient k, is actually fairly consistent with experimental data for kinesin [T] I32| . 
provided the tail is not stretched too far ( r < 70 nm) |21| . and we use it for clarity in our main development. 
We stress that we consider the linear force law ^ merely as a convenient and reasonable phenomenological 
approximation; the connecting tail has an effectively jointed structure which gives rise to a more complex 
relationship between the force and the motor-cargo separation. Our analysis can be readily extended to 
nonlinear spring models that have nonzero rest length [57], a sigmoid function as in [2S] and the wormlike 
chain model | i25) , provided the force law is approximately linear over a wide enough range ( [A| . 

2.3 Force-dependent instantaneous velocity 

We describe next the nondimensional drift functions g, which summarizes the detailed response of the motor's 
transport properties to an applied force. The function g can be thought of as a nondimensionalized version 
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Motor Properties 



Quantity Label Value for kinesin in vitro 



l^^^mm^m step size 


L 


nm ^mm^^^^^^^i 


Stall force 


F, 


7 pN 


Unperturbed motor velocity 


V 


nm/s ^^^^^^1 


Neck linkage spring constant 


K 


0.34 pN/nm 


"Randomness Parameter" 


r = a'^ /Lv 




Effective motor diffusion 




5000 nmVs 


Cargo and Environment Properties 


Effective Cargo Radius 


a 


500 nni 


Fluid Viscosity 


V 


1 X 10"^ pN • s/nm^ ^ 


Optical trapping force 


e 


-10 to 20 pN 


^Boltzmann constant • temp 


ksT 


4.1pN-nm '^^^^^J 


Cargo friction 


7 — Givar] 


1 X 10"'' pN ■ s/nm 



Table 1: Experimental estimates of physical parameters for kinesin |671 [5U1 132|. Our estimate for the 
cargo friction is based on the viscosity of water. In Section |4] we assume that the in vivo viscosity can be as 
high as 6 X 10~^ pN • s/nm^. consistent with |43| . 



of the force- velocity curves typically reported in experimental measurements of motors |^ EH E2 ES] , with 
the unforced motor velocity and stall force scaled so that g{0) — 1 and g{l) — 0. Note that the argument / 
of g{f) is the ratio of the applied force, measured in the direction opposing natural motion of the motor, to 
the stall force. In general, these functions g decrease at an order unity rate for substall forces < / < 1, and 
they saturate at finite constants both for superstall forces (/ > 1) and forces applied along the direction of 
the motor's natural motion (/ < 0). Theoretical models |46I[M| generally employ simple functions satisfying 
these conditions, though some models actually incorporate, for tractability rather than accuracy, functions 
g that grow linearly for large applied forces / |46| . As a specific example that we employ in numerical 
simulations, we take a sigmoid function form for g, namely 

g{f) = A-Bt^iMCf-D). (3) 

We impose the constraints ^(0) = 1 and ^(1) = 0, as well as the asymptotic relations: 

hm g(f)^ , hm 5 / = , 

with Wmax the maximum speed the motor moves when pulled by a strong assisting force; Wmin the nonpositive 
velocity at which the motor moves backwards when dragged by a strong superstall force [6J; and v the 
unperturbed motor velocity. This uniquely determines the constants to be 



* ^max ^~ ^min 7-, ^max ^min , i ' 'i^ma.x 

A = , B = , D = arctanh 

2v 2v 



and C — arctanh — arctanh 



2v 



Such a function with physically relevant parameters is depicted graphically in Figure [2j Although we have 
described a specific model for the force- velocity relationship, our methodology and conclusions apply to 
more general choices of models. As we assert in Theorems |5.1| and |5.4[ the concavity of g near zero force 
and near stall force is critical to qualitative behavior of a motor-cargo system when more than one motor 
is attached to the same cargo. (This is also noted by Wang and Li [72] and Kunwar and Mogilner |36j). 
Our operating assumption is that the force-velocity curve is concave down for small external forces, and 
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Example of a force -velocity relationship. 




Figure 2: An example model ([3| for the force-velocity relation vg[f / F^) with maximum velocity Wmax = 600 
nm/s, minimum velocity fmin = ^50 nm/s, free velocity v = 500 nm/s and stall force F.^, = 7pN . 

concave up for near-stall forces. This corresponds to typical force- velocity curves for molecular motors under 
normal operating conditions |67l [B] , but we note that the force- velocity data from these sources also indicate 
a concave-up behavior for small external forces when the local ATP concentration is very low. (See H5] 
for further discussion and contrasting properties for dynein). We summarize these qualitative properties, 
typical of the high ATP regime, as follows. 

Assumption 1. We assume that the instantaneous force-velocity relationship g : M — > M is smooth, and 
hounded |^ with hounded derivatives. We normalize g with respect to the unencumhered velocity of the motor 
and the stall force such that 

5(0) = 1, andg{l)=0. 
Furthermore we require g to satisfy the following qualitative properties: 
(i) (Monotonicity) g{f) is strictly decreasing in f f67}j : 
(a) (Concavity) There exists an f^ G (0,1/2) such 

9"{f ) < 0, for all f e (-oo, and 
g"{f)>0,for a///e [l-/,,oo); 

(Hi) (Strong Concavity) The function ri{f) :— g{f) + g{~f) is decreasing in f for f > 0, while f}[f) := 
g(l + /) + g{l — f) is increasing in positive f. 

Note that functions of the form ^ satisfy Assumption [l] Indeed, g'{f) = —BC/ cosh^(C/ — D) and so 
when |wi„in| < l^maxl and v > (wmax + Wmin)/2 (so that A, B,C,D > 0), 

(<?(/) + 9{-f)) = -BC ( I \ \ < 

q/ \cosh {Cf-D) cosh {-Cf-D) J 

implying part the strong concavity property about 0. A similar calculation confirms the property with 
respect to / = 1. 

We note though that the force-velocity relationship in the data collected by Carter and Cross [6, is not 
monotone for negative values of /. In fact, it is difficult to detect the precise properties in this regime because 
the motor is so apt to unbind from the microtubule. We choose to take from the data that the velocity of 
the motor is bounded, which is consistent with the observation that the stepping must be rate-limited by 
the binding rate of the forward head of the motor. 
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3 Nondimensionalization and Identification of Distinguished Limits 
for Asymptotics 

We begin with a consideration of the low viscosity environment characteristic of in vitro experiments, for 
which the physical parameters in our model take values approximated by those presented in Table |2.2[ The 
equations ([l]) can be readily simulated, but we exploit asymptotic techniques to reduce the coupled dynamical 
system to a more easily interpretable lower-dimensional system. To perform a nondimensionalization, we 
first choose suitable reference length and time scales to render the system ([T|) nondimensional and to identify 
small parameters. For this approach to produce physically meaningful small nondimensional parameters, the 
reference units must be chosen so that the nondimensional variables vary on order unity scales |401l50]. For a 
cooperative system of motors, the length and time scales can be chosen in the same way as for a single motor 
attached to the cargo, provided is not small compared to the resulting small nondimensional parameter. 
This latter consideration is well-satisfied, since < 10 and the small parameter that will emerge from the 
analysis has magnitude e ^ 10~'^. We remark that in other configurations, care is needed. In particular, 
in a tug-of-war configuration, the motors are all nearly stalled, and this will induce different length and 
time scales to their dynamics as compared to the fully cooperative configuration. Moreover, while our 
nondimensionalization is still valid for the higher viscosities characteristic of in vivo cellular environments, 
our conclusions about the relative sizes of the nondimensional parameters (particularly e) will need to be 
revisited for this setting. 

We consider the relative size of the forces acting on the motor-cargo system, using the physical magni- 
tudes of the various parameters as listed in Table [23] characteristic of an in vitro experiment. The frictional 
force 7u that would be exerted on the cargo if it were pulled at the full unencumbered speed t; of a mo- 
tor is 2.5 X 10"'^ pN, by far the smallest force incorporated in the model. To characterize the strength of 
thermal effects, we suppose that the tail connecting the motor and cargo is in approximate thermal equi- 
librium, which implies (under the linear spring model) that its distortions impart a typical potential energy 
^fceT, so that ^kE [(X(i) — Z{t))'^] ~ ^kuT. Consequently, the length of the spring has mean-square value 
E [X{t) — Z(t))^] ~ k^T/K, so the typical magnitude of force fluctuations induced for the connecting tail is 

,1/2 



kE [{X{t) - Z{t)f] ' - ^Jk^Tn - 1.3pN. 

The stall force required to arrest the motor is, depending on the actual motor, in the range of 5 to 10 pN. 
In experiments, the laser traps are generally used to apply forces that range from —10 to 20 pN. 

These suggest the following leading order picture: In the absence of a laser trap, the forces induced by 
thermal fluctuations are substantially stronger than those induced by friction, so the dynamics of the tail 
connecting the motor and cargo should be dominated by thermal fluctuations. We therefore infer that a 
suitable reference length scale is ^Jk^T / n ~ 3 nm, the typical length of the fluctuations in the tail due to 
thermal effects. The connecting tail itself, in a realistic model, is actually on the order of 70 nm, but as 
we discuss in[^ this length scale does not play any dynamical role, so we do not use it as a basis for our 
dimensional analysis. Also, the 8 nm step size of the molecular motor may appear relevant, but we are 
considering the effective behavior of the motor-cargo complex as it progresses over much larger distances; in 
our model this discrete-stepping detail has already been coarse-grained, with the effects incorporated into 
the force- velocity and force-diffusivity relations from Subsection |2.3[ 

To estimate the time scale, we note that neither the thermal forces nor the friction force are particularly 
strong compared to the stall force, so we crudely approximate the force-velocity relationship by its local 
linearization about zero force: 

g(/)«l+5'(0)/, 

where we have noted that f;(0) = 1. Under this rough approximation, the dynamics of a cargo with a single 
motor can be written as the linear system: 



X(t) 

m 



= J 



X[t) 

m 



dt 



-Oh 



dt 







^2kBTh 



dW^ii) 

dwM 



(4) 
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with constant matrix 



J 



(5) 



F, y F, 

K K 

7 7 

This matrix has eigenvalues (corresponding to no resistance when the motor and cargo move together) and 



A = -^^-^, (6) 

7 

whose magnitude describes the effective rate constant of the tail stretch X{t) — Z{t). The inverse of its 
magnitude provides a good reference time scale. As noted above ^ ^v, so A — ^ and we choose ^/ k as 
the reference time scale. This is precisely the time scale which characterizes the dynamics of the tail when 
attached to the cargo and a fixed object, which is appropriate, since the motor dynamics are relatively slow 
compared to the cargo, a fact we will see in the sequel. This choice of reference time scale is also valid for a 
wide range of laser trap forces and for the case of multiple cooperative motors. It is true that for laser trap 
forces large enough to bring the motors near to a stall, the length and time scales change, but our aim is 
to select reference length and time scales that work well for a range of laser trap forces, from zero to stall. 
By neglecting the laser trap force, we can specify length and time scales that are relevant for a much wider 
range of forces than those based on stalled configuration dynamics. 

While mass units technically enter an equation involving force balance, none of the variables of interest 
involve mass units and we can simply divide the equation for Z{t) in Eq. [T| by 7 so that mass units dis- 
appear from the left hand side and therefore also from every group of parameters on the right hand side. 
Consequently, we effect our nondimensionalization by combining this simple manipulation with the explicit 
rescaling of variables with respect to length scale yj2k^T/ k, (the factor of 2 for convenience) and time scale 

xii) = zii) ^ 

where t — {ti/^)t. We also employ an important fact about Brownian motion: that scaling time, W{at), is 
equivalent in distribution to scaling space, ^/aW{t), for any positive a. The governing equations in terms of 
these nondimensional variables therefore become 

dX,{i) = eg{s{Mi) - Z{i))) dt + (T,„/e dWS) (7) 

N / n\ 



1=1 



where we have identified the nondimensional parameters listed in Table |2] 

We find that, in vitro, the smallest of the nondimensional parameters is e '-^ 3 x 10~^. The diffusion 
coefficient ratio a^/^ ~ 8 x 10~^ is also small relative to the parameters s and 9. This suggests the bio- 
physical relevance of the following distinguished asymptotic limit, which we will see also has mathematically 
desirable properties: e is taken as the fundamental small nondimensional parameter, and we rewrite (Tm/c 
as a parameter proportional to ^/e, i.e., am/c = ^/eP' '^here p is defined in Table |2] The nondimensional 
parameters s, 9, and p will be treated as order unity and independent of the small parameter e. 

As we have noted, the viscosity of the fluid environment ranges widely from the experimental setting to 
what is expected in the cytoplasm We therefore treat 7 as a flexible parameter, and our distinguished 

asymptotic limit described above remains relevant as 7 is varied (noting how it appears in the nondimensional 
groups in Table[2]), provided the parameter e remains small, which we shall see in Section 4.1 means no greater 



than 0(10""'^). If so, we call the dynamics diffusion dominated: otherwise for reasons that become clear in 
Section [4] we call the dynamics drag dominated. 

With the distinguished limit described above, the small parameter e appears in the governing equations 
in two places in the equation for the motor variables Xi(t). Because of the appearance of e in the drift 
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Dimensionless Groups 



Description 


Label 


Definition 


Value in vitro 


Ratio of the friction force exerted by the 
cargo when pulled by the natural motor 
speed to the typical force exerted by the 
tail due to thermal fluctuations. 


e 


vy 
y/2kBTn 


3 X 10"^ 


Stallibility. The ratio of the typical tail 
force to the stall force. 


s 


y/2kBTK. 
F, 


0.2 


The ratio of the force applied by the laser 
trap to the typical tail force. 


e 


6 

y/2kBTK. 


—6 to 12 


The ratio of the effective diffusion coefficient 
of the motors to the standard Stokes- 
Einstein diffusion coefficient for the cargo. 


2 


2kBTh 


6 X 10"^ 


Randomness parameter 68i 1111 
or inverse Peclet number \35\ 
with respect to length scale of 
thermal fluctuations in tail. 


P 


v^2kBT 


2 



Table 2: Nondimensional groups with the range of values associated to the constants in Table |2.2| with 
viscosity corresponding to that of water (representing an in vitro environment). 



term and ^/e in the diffusion term, we can interpet as a "slow" nondimensional time scale characterizing 
the motor dynamics relative to the order 1 "fast" time scale characterizing the cargo (or equivalently, the 
connecting tail). The physical justification for this time scale separation is that the cargo friction is very small 
relative to the spring force, and therefore the cargo (and spring) equilibrate quickly relative to the motor's 
motion along the track. The relevance of thermal fluctuations, however, implies that this equilibration is 
not to a rest point, but rather to a stationary probability distribution of fluctuations about the minimal 
energy point of the spring. More speciflcally, consider the time change t = €~H. Application of stochastic 
averaging theory (see [T2l 165], |29] ) allows the derivation of effective equations for the motor coordinates, 
without explicit reference to the cargo position (see [541 1551 [TU] for analogous reductions for single-motor 
case). When e is small, the behavior of the motor dynamics {Xi(e~^f)}"^j^ is well-approximated, in the sense 
of weak convergence, by the solutions {Xi{t)}'^^^ of the following SDE, 



where 



■jTz{z;yi,e) = -\/ — exp 



g {s{xi - z)) 7rz(x, z; 0) dz; 



'N z 



Y.1 _ ^ 

N N 



(8) 

(9) 
(10) 



The function ttz is the density of the stationary distribution of the cargo position, given fixed motor positions 
X = {xi, . . .,xn}- 

Note that while the cargo variable has been removed from explicit consideration, the effective dynamics 
of the motors are coupled together because they are all connected to the cargo. That is, the effective drift 
coefficient gi for motor i depends on the current positions of not only motor i but the other motors as well. 
By following the changes of variables we have performed in this section, we conclude that the simplifled 
equation (|8| provides a good approximation (provided the parameter e is sufficiently small) to the original 
model equations ([T|) in the sense that the statistical dynamics of the solutions {Xi{t)}fLi to the true model 
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equations ([T|) are well approximated by the statistical dynamics of yj2k^T / kX{J^ evaluated at the rescaled 
nondimensional time t — vyj nj {2kBT)t, where {Xi{t)}fLi are solutions to the simplified equations ([8|. We 
will explicitly check the quality of this approximation in Section |4] 



4 Analysis of transport by one motor. 

Numerous authors have expressed concern about translating experimental observation into in vivo predic- 
tions. Shubeita et al [64 and Mitchell & Lee |32| provide excellent discussions on the topic. At issue is 
that the now classical force- velocity relationships [67! ^ were measured in the presence of a water-like fluid 
environment, while the viscosity in vivo may be 100-1000 times greater. To infer the effect of high viscos- 
ity, Shubeita et al [64] attempted to measure the stall force by trapping lipid droplets that are driven by 
kinesin-1 motors. While they measured the low viscosity stall force to be near 5 pN, the apparent stall force 
in vivo was closer to 2.4 pN. In a more theoretical treatment, Mitchell and Lee |33] used the Stokes-Einstein 
relationship to estimate the drag force incurred on a motor via various shapes and types of cargo. Their 
goal was to estimate whether a single motor could overcome the anticipated drag force in order to produce 
experimentally observed motor speed. 

As was indicated in the previous section, we flnd that there is a qualitative change between the low 
and high viscosity regimes. We will show that nevertheless our modeling framework can provide a smooth 
interpolation between them to explain how the presence and dynamics of the cargo affects the speed and 
effective stall force for a single motor. 



4.1 Low viscosity regime and inference for g 

When the viscosity of the fluid environment is similar to that of water, we can compute the mean velocity 
and diffusivity of the system by averaging the force- velocity curve g against the quasi-stationary distribution 
of the cargo, ttz, with the motor position essentially held constant. For the linear spring model we are using 



in the main text, the stationary distribution of the cargo is Gaussian (10 1, and the average velocity and 
effective diffusivity for one motor with cargo can be expressed as 



t^OO t 



Vi{e) = lim - 

t—^oo zt 

where we have introduced the Gaussian average: 



(11) 



g{s{x ~ z))^e 



dz = A[g] (sO, y ) . 



(12) 



The last equality follows by introducing the variable y = s{x — z). As expected, this effective drift Vi{9) is 
not dependent on the location of the motor, only its averaged distance from the cargo. 
Written in terms of the original parameters, we have 



(13) 



Since ^^i^ = 0.03 <C 1, it follows that for a given external force 9, Vi{6) w vg{9/ F^), which corresponds 
to the naive estimate that the motor should move overall at a speed corresponding to the applied force 9 
and the force-velocity relationship g. 

Furthermore, this characterization of the average velocity gives a clean representation of the effects of 
perturbing the physical parameters of the system. For example, the velocity given by ( 13 1 does not depend 
on the viscosity of the fluid environment, consistent with observations in vitro [16 and in silico [36'. The 
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stability of the velocity with respect to viscosity encourages the notion the observed force-velocity curve 
should be robust as long the assumption that the ATP concentration remains constant holds. 



The quantity (13 1 is the experimentally observed velocity, so the theoretical instantaneous force- velocity 



curve g is the function that best fits Vi{9) after "inverting" the Gaussian average. Such inference can be 
done using the recently developed non-parametric statistical methods of Asencio et al [T]. Figure |3] compares 
the instantaneous force-velocity curve for a motor with the "dressed" version that accounts for the effects of 
thermal fluctuations in the cargo (and therefore connecting tail). 



A second example of a prediction made by ( 13 1 is 



found by perturbing the spring constant k. For small 
9, we note that the velocity decreases when the spring 
constant increases. This is due to the downward con- 
cavity near zero and the fact that increasing k widens 
the Gaussian average. 



Effective Force-Velocity Relation 
N=1 Motor 



One final novelty of the formula ( 13 1 is that if g 




-5 5 

Laser trap force (pN) 
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were concave up near 0, as is the case in the low ATP 
regime, the opposite inequahty would hold, meaning 
that the fluctuations would improve motor progress. 
This can be interpreted as an example of the phe- 
nomenon called strain gating [37]. It has been ob- 
served that when force is applied to the motors, they 
can take steps even when little or no ATP is avail- 
able. In the presence of cargo fluctuations, the cargo 
will explore regions of space that aid the motor in tak- 
ing steps. Because the motor would respond more to 
pulls from the cargo when it fluctuates in the forward 
rather than the backward direction, the average ve- 
locity of the system is greater in the forward direction 
than it would be in absence of cargo fluctuations. 

All of these properties hold as long as the time 
scale separation holds between the motor and cargo 
dynamics, i.e., as long as e is relatively small. In Fig- 
ure [3] we present the average velocity of the original system as a function of the cargo friction, which is 
proportional to the solvent viscosity. In agreement with a similar plot in |36| . the velocity holds approxi- 
mately constant for two orders of magnitude greater than that of water, or in our nondimensional terms, for 
e < 0.1. Beyond that we must use a different analysis, which is provided in the next section. 



Figure 3: Comparison of the effective force- velocity 
relationship of a single motor attached to a cargo 
from direct stochastic simulations of the model equa- 
tions ([T| (blue with error bars denoting sample stan- 
dard deviation) and the stochastic averaging result 
(12 1 (red solid). We also show (magenta dashed) 



the instantaneous force-velocity curve for a motor 
vg{6/F^) which neglects thermal fluctuations. The 
curves are somewhat difficult to distinguish as they 
are almost coincident. Parameter values are listed 
in Table [221 



4.2 Crossover to high viscosity regime 

Because the stochastic averaging no longer holds when the solvent viscosity (or bead size) is so large that 
e ^ 1 or greater, we return to the system (|7|. In the one motor case, there is a simplification that arises 
from introducing Y — X ~ Z , the time-dependent distance along the microtubule between the motor and 
cargo. After a nondimensionalization and under the timescale t = e^^t, we obtain 

dX{t) = g{sY{F)) dt + ^pdW^it) 

dYit) = [g{sYit)) - -^{Yit) - e)] dt+ ^/^TlTedWyit). (14) 

where Wy = ^ {y^Wx — \f\J(Wz) is a standard Brownian motion (correlated with Wx\ We see that 

the drift term of y is a function of Y alone, and the diffusion term is simply a scaled Brownian motion. Now, 
as a one-dimensional SDK, the drift term for Y can be written as the derivative of a potential function, and 
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we can solve for the stationary distribution of y as an exponential of this potential: 



T^viy) = C exp 



1 



-{y~fff + 2e / g{sy')dy' 



(15) 



This formula reveals the manner in which the high viscosity regime is a perturbation of the low viscosity 
regime: when e is small, Try is essentially a Gaussian density with mean 9. As before, to find the long-term 
velocity we average the drift term for X against this stationary density. 



hm = hm = v / g{sy)nY{y)dy. 

t-i-oo t t-i-oo t 



(16) 



This expression can be rewritten more suggestively as: 



lim = V 

t— >oo t 



l + ep 
2e 



dTTviy) , 2iy-e) 



dy 



l + ep 



TTviy) 



dy 



{y - e)TTY{y)dy 



v{{Y) - 9) n{Y) 



7 



(17) 



where (•) denotes a statistical average over the fiuctuations in the connecting tail length (governed by 7ry(y)). 
The final expression indicates that the speed of the system is given by the average net force on the cargo 
(connecting tail force minus applied load) divided by its friction coefficient. Note that the nondimensional 
expressions approach 0/0 in the e -> limit, which is why we did not recast our stochastic averaging 
expression ( 13 1 in this form. 



To gain some further intuition, it is worth considering a deterministic version of the system (14) defined 
by the system of ODEs, 



di 



To find the average velocity, we find the stable equilibrium for Y{t), namely a value y^ such that g{sy^,) = 
j{y* — ^)- Because g is decreasing in y*, the right-hand side is of g^^"*-''^ is strictly decreasing, so we know 
that if an equilibrium exists, this point is unique. For the moment, assume that g is linear for y S [0, 1] and 
therefore of the form g{y) = 1 — y to maintain consistency with Assumption [T| Then the above equation 
becomes 1 — sy* = (y* — 0)/e and solving for y.^, gives y^, = (e + 9)/{es + 1). 

This implies that, for 9 € [0, 1], the asymptotic nondimensional velocity for linear g satisfies 



^-X(^)-5(.sy,) = ^(y, 



1 - s9 
1 + se ■ 



Note the similarity to the penultimate expression in Eq. (17 1, with the statistical average over fiuctuations 



now replaced by a deterministic value for the stable length of the connecting tail. In the original parameters 
of the problem, we have the approximation 



Vi{9) 



v{l~9/F, 
1 + jv/F, 



for 6'e [-F,,F^]. 



(18) 



In fact, this approximate formula remains unchanged even in the presence of thermal fluctuations (p > 0), 
as can be shown by by simply substituting g(y) = 1 — y into (16 1. Figure |4] shows that this simple explicit 
formula accurately represents the effective behavior of a single motor even into the high viscosity regime 



(beyond the validity of the stochastic averaging formula (11 1), for applied forces up to 5 pN in magnitude. 

In particular, this formula gives a simple representation for how the velocity of a single motor with zero 
external force should decrease as a function of viscosity. We immediately read off a critical nondimensional 
ratio, jv/F^,, which compares the motor's stall force to the drag force that would be felt by the cargo if the 
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Unloaded Velocity 
N=1 Motor 

550 1 ' ' 




Cargo friction y(pN s/nm) 



Effective Force-Velocity Relation 
N=1 Motor 

800 1 ' ' ' 1 




Laser trap force (pN) 



Figure 4: Performance of the linear approximation. Left panel: Velocity of cargo attached to single motor 
as function of cargo friction parameter 7. Results of direct simulation of ([T]) are indicated by the blue curve, 
with error bars representing one sample standard deviation. The black dotted curve represents the linear 
approximation formula (18 1, whereas the red dashed line indicates the theoretical value corresponding to 
the e 4, limit. The values of 7 presented correspond to a bead of size a — 500nm [TF and solvent viscosity 
ranging from the viscosity of water (at 300 K) to 1000 times that value. Right panel: Comparison of the 
effective force- velocity relationship of a single motor attached to a cargo with friction 7 = lO^^pN • s/pN 
(e = 3) from direct stochastic simulations of the model equations (jlj (blue with error bars denoting sample 
standard deviation), the precise theoretical formula (16 1 with the high viscosity modification (red solid), 
and the explicit formula ( |18[ ) for high viscosity modification based on a linear approximation (black dotted). 
We also show (magenta dashed) the instantaneous force- velocity curve for a motor vg{6/F^,) which neglects 
thermal fluctuations. Other than the cargo friction 7 (and proportionally related viscosity 77), parameter 
values are listed in Table [2T2l 



motor could achieve its natural low viscosity speed. As can be noted from Figure |4| the changeover from 
this ratio being negligible to significant occurs near the in vivo viscosity estimated by Mitchell & Lee [43j . 

The formula ( 18 ) does not do well, however, at large applied force; it predicts the stall force to remain 
exactly F^,, regardless of the viscosity of the fluid. As we will show in Section 5.2 the stall force of the full 



stochastic model in (16 1 remarkably increases with the viscosity of the fluid due to the upward concavity of 
g near stall force. This claim is counter to the observations of Shubeita et al [M], for example. A deeper 
version of the model and further analysis are necessary for making predictions in the near-stall regime. 



5 Analysis of transport by multiple identical motors 

Experimentally it is very difficult to determine precisely how many motors are attached to an observed cargo 
at a given time. Typically, in order to get some certainty, transport is observed in a fluid medium that has a 
very low concentration of kinesin. A majority of the observed cargos are immobile (or are undergoing simple 
diffusion) but a minority are seen to move almost deterministically for a time. Because the kinesin are so 
rare in the solution, it is hypothesized that only a single kinesin is bound to the cargo and a microtubule. 
When the concentration of motor proteins is increased, not only is the mobile portion of the cargo population 
increased, but the run length and response to external forcing is changed. In particular, motor concentration 
is expected to be often large enough in vivo for multiple motors to bind to a cargo. This motivates the 
development of quantitative relationships between in vitro observations involving a single motor bound to 
the cargo to in vivo predictions in which multi-motor transport is presumed to be relevant. The purpose 
of most of this section is to derive such analytical connections within the modeling framework laid out in 
Section [2 

A striking observation from theorists and experimentalists is that the attachment of multiple motors to 
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a given cargo does not necessarily imply a larger average velocity of the motor-cargo system. This has been 
observed separately by Kunwar et al fST, Wang and Li fT^I and Bouzat and Falo Interestingly, each 
group gives a distinct reason for the phenomenon: Kunwar et al cite detachment and attachment dynamics; 
Bouzat and Falo cite excluded volume effects among the motors; and Wang and Li cite the structure of the 
force-velocity curve. We highlight the multiplicity of explanations in order to emphasize the richness of the 
underlying biophysics and the benefit of careful analysis. Our model also reproduces the two-slower-than-one 
result in the low viscosity environment, for reasons closest to that of Wang and Li. However, as we see in 



Section 5.4 our model predicts that this phenomenon disappears when viscosity is closer to that expected 
in the cytoplasm. 



5.1 Analysis for a semi-deterministic approximation 

The two sources of stochasticity, cargo fluctuations and variation in the spatial configuration of the motors, 
both contribute in a non-trivial way to two-motor transport. When both are eliminated, we have the model 
for force distribution introduced by Miiller, Klumpp and Lipowsky |45 l l45 l l47 ) . Each motor is presumed 
to bear half the load of the cargo and the resultant drift term in the original parameters of the system is 
vg{6/2F^). We will refer to this as the force-balance theory. 

In order to distinguish between the contributions of the two sources of randomness, we add one back 
and then the other. Because it turns out to account for a dominant contribution of new behavior, we first 
analyze the system with changing spatial configurations in the absence of cargo fiuctuations. Versions of 
such an approximation, which is not a physically natural asymptotic limit of the system, have been explored 
before [721 E5] . Our analysis differs from prior work in our SDE approach to describing the evolution of the 
configuration of the motors over time. 



(xi, X2, ■ ■ ■ , Xn), its mean position is Z[x) = 



By inspection of the cargo's stationary distribution ttz from (10 I, for a given motor configuration x = 



For our simplified model with non-fiuctuating 



cargo, we simply assume the cargo instantaneously adjusts to this mean position as the motor coordinates 
evolve. We denote the dynamics of the motors with non-fiuctuating cargo by {-^(f)}ie{i,Ar}, which satisfy 
the system of SDE 

dl,(t) = g{s[X,{t} - Z{±m)dt + VpdW.it). (19) 

For a two motor system, the cargo position is given by Z{x) — {xi + x-2)j2 — 0/2. It follows that the 
nondimensionalized opposing force felt by the first motor is s{xi—Z{'k)) — s{xi—X2+0) /2. The opposing force 
felt by the other motor is s{x2 — Z{'x)) = s{x2 — xi + 9) /2. To simplify notation, we introduce a new function 
G(r) :— g{—sr/2) which represents - in the zero external force setting - the instantaneous drift induced on 
a motor when the other motor is a given signed distance away. That is, letting r — xi — X2, (which is to say 
that X2 is a signed distance — r from xi), the drift experienced by the first motor is g{s{xi — X2)/2) =^ G{-~r). 
When there is an external applied force 9, the drift of the first motor is g{s{xi — X2 + 6)/2) = G{—r — 0) 
and the drift of the other motor is g{s{x2 — xi + 0)/2) = G(r — 9). 

By changing variables to center-of-mass M{t) — }^{Xi{t) -\- X2{t)) and difference R{t) — Xi{t) — X2{t) 
coordinates, we obtain 



dAf(f) = - G(R{T) -§) + G{-R{t) 



dR{t) = - G{R{t) -9)- G{-R{t) - 6) 




2pAWr{t). 



(20) 



where we have introduced Wm = '^(Wi + ^2) and Wr = ^(W^i ~ ^2), which can be treated in law as 
independent standard Brownian motions. 

The process R is independent of M and the drift term can be viewed as the derivative of a potential 
function. Therefore, introducing the notation G±{r; 9) = G{r — 9) ± G{—r — 0), we can express the density 
of the stationary distribution of motor separation process R as 

1 r 1 

/ G^{r';9)dr' , -cxX r < 00, (21) 
P J -00 
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where Cr is a normalizing constant. The midpoint process M{t), which we use as a proxy for the position 
of the entire system, is then a functional of the separation process R(t) plus a Brownian motion, 



Applying the law of large numbers and the ergodic theorem for diffusion processes, we can calculate the 
asymptotic average velocity limf_j.oo M{t)/t by integrating the drift of M against the stationary distribution 



M{t) 



lim 

t^oo t 



G+{r;e)7rRir;0)dr. 



The asymptotic diffusivity is given as follows 



— OG \J —CO 



— dr. 



(22) 



(23) 



This effective diffusivity is obtained by way of the central limit theorem for stationary stochastic differ- 
ential equations We center the process M about its mean and rescale by the square root of time, and 
use Eq. (5.11 to write: 



^_ (a/ (t) - v,{0)t) = i,^* ^G+(i?(i'); 0)) ~ v,{e)dt' + 



(24) 



The first term on the right converges in distribution to a Gaussian random variable with mean zero and 
variance 



{\G+{r'-e)~V2{d)>R(r'-,0W 



1 



— dr 



by the central limit theorem given in Chapter 1 of [38J. 

Together with the Brownian motion term in (24 1, which is independent of R{t), we conclude that the 
external force dependent diffusivity is given by ( 23 ) . Thus for a given large i, 



where Z is a standard normal random variable. 



2'b2{e)tz, 



5.2 Comparison of One-Motor and Two-Motor Systems 

We are now ready to characterize how the spatial distribution of the motors affects the effective transport 
properties relative to the simpler models |151 H5] in which all bound motors are assumed to share the load 
equally. 

When there is no external forcing, we can understand the interference between the motors via the following 
thought experiment: Suppose the motors are separated by a distance r and that the cargo is fixed at the 
mid-point between the motors. Each spring connecting the motors will be stretched equally inducing an 
opposing force with signed magnitude of ±sr/2 on the leading and trailing motors respectively (positive 
opposing force on the leading motor; negative opposing force on the trailing motor). The force- velocity 
curves in the literature |^ |321 E21 ES] seem to indicate that a "helpful" force can only speed up the motor a 
small amount, whereas the same force applied in the opposite direction can slow the motor down considerably. 
Therefore the mean of the two induced velocities is less than the velocity of a single motor, as characterized 
by the inequality 

^(g(rs/2)+.g(-rs/2)) <.g(0), r ^ 0. 
This is the motivation for emphasizing the concavity properties of g in Assumption [T] 
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In the presence of an external force 9, the average velocity of the two motors will be less that that of one 
motor when 9 and the separation distance r satisfy 

^{g{s{r + §)/2)+g{s{-r + 9)/2)<g{s9). (25) 



For the class of g functions typical of force- velocity curves (such as our schematic example in Subsection [2^ , 
this inequality will only hold for relatively small values of external force 9 and a possibly restricted range of 
separation distances r. For example, note that when r — 0, the left-hand side is g{s9/2) which is actually 
greater than g{s9) when 9 > 0. Nevertheless, we must average the left hand side of ([25| against the stationary 



distribution of the motor separation r (21 1, to attain a comparison of the average velocities of one-motor 
and two motor systems. This explains the restrictions in conditions under which we can be confident that 
two motors bound to a cargo will be slower than one. 

Proposition 5.1 (Two motors can be slower than one). Suppose that the force-velocity function g satisfies 
Assumption [7} Then for any given value of the ratio of typical spring force to stall force, s, there exists a 
positive value 9c{s) such that for all 9 < 9c{s), 

V2{9) <Vi{9). (26) 



Proof. In light of ( 13 1 , the nondimensionalized velocity of a system with one motor is given by Vi {9) — g{s9) . 



By ( 22 1 and the definition of G, the inequality ( 26 1 holds whenever 

1 



2 ,g{s{r + 9)/2) + g{s{-r + 9)/2))dTTnir; 9) < g{s9). (27) 



As discussed earlier, the inequality (251 does not hold for all r, so we must establish the inequality for 9 = 
and then rely on continuity of the functional with respect to 9. Indeed, we employ Assumption [l| ([iii|) to find 
that 

1 

-{g{sr/2) + g{-sr/2))d7rR{r; 0) < 5(0)d^fl(r; 0) = g{0). 

To extend the inequality to an open neighborhood of the origin, we note that g{s9) — V2(^) is continuous in 
9 and so must remain positive on some open interval containing the origin. □ 

At large applied forces, the appropriate comparison is between the stall force of two motors versus twice 
the stall force of one motor. First it is worth recalling that in nondimensional terms, the stall force is 
9 = 1/s (see discussion in Section [s]). Therefore we compute the velocity for 9 = 2/s and demonstrate that 
it is positive. Indeed at a separation r, the velocity of the system satisfies 

^{g{s{r + 9)/2) + g{s{-r + 9)/2)) = i(.g(l + sr/2) + g{l - sr/2)) > g{l) 

where g{l) = by hypothesis. 

The intuition here is that while the leading motor is actually moving backward because it is experiencing 
a force greater than its stall force, the trailing motor is still moving forward, and doing so with a magnitude 
greater than that of the leading motor. The mean of the two velocities is therefore positive and the system 
has an overall positive drift. In other words, it takes more than twice the stall force of a single motor to stall 
the two-motor system. 

Proposition 5.2 (Superadditive stall forces). Suppose that the force-velocity function g satisfies Assumption 
[7] and that the stallibility satisfies s < 1. Let 91 and §'2 denote the minimum nondimensionalized force 
necessary to result in an average velocity of for a cargo attached to one or two motors, respectively: 

9* = mm{9 : Vj{9) < 0}. 

Then 

91 > 291 (28) 
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Proof. For a given 9, we recall Vi{0) — g{s9). By Assumption 1, we can directly observe that 61 = 1/s. To 



show (28) recall the definition of 9i from Assumption |l | (|iii[) . It follows that for all 9 > 29i/s 



W)^ I \{9{s{r + 9)/2)+g{s{~r + 9)/2))d7Tn{r;9) 



> / g{.s9/2)dTrRir;9)=gis9/2) = V,{9/2) 

JR 

In particular, since s < 1 we have 9i < 1 < l/s — 91. Therefore for all 9 E {29i/s,9l], we have V2{0) > 



Vi (61) > implying (281. □ 



Because the the strict superadditive inequality relies on averaging against the stationary distribution ttr 
of the motor configurations, we see that stochasticity in motor position is crucial to this result. Without 
these fiuctuations, the motors would simply stall at 9 = 2/s because each motors' share of the total force 
would be exactly their stall force. The random component of the stepping makes it possible for the motors 
to get out of a stalled configuration. Once separated, the differentiated responses of the leading and trailing 
motors makes progress possible for the motor-cargo system. 

The condition s < 1 means that the typical spring force is less than the stall force of the motor, so that 
the fluctuations in the cargo by themselves do not typically subject the motors to stalling forces. We remark 
that similar conclusions relating the properties of the two motor and cargo system to the convexity properties 
of the single-motor force-velocity curve, also for the case of non-fluctuating cargo, were developed by |72| in 
the context of an exact solution of a model in which the motors were represented as random walkers on a 
lattice, coupled by their attachment to the cargo. Their model did not include cargo fluctuations, but in the 
next section we will extend our analysis to include these. 

Before moving on, we note a contrast to models that do incorporate motor unbinding. As discussed 
earlier, part of our objective is to determine which properties of multi-motor transport can be directly 
attributed to distribution of forces when the spatial configuration of the motors is accounted for. It has been 
shown, however, that in certain regimes motor detachment is the dominant feature. 

For example, Kunwar et al [37, 36, observed that when the cargo is subjected to a significant force, the 
majority of this force is borne by the leading motor, which is therefore more likely to detach than the trailing 
motor. When detachment occurs, the system snaps back to the location of the trailing motor. Often this 
other motor, now bearing all the force of the optical trap, will detach as well. (This detachment cascade is 
crucial to the tug-of-war dynamics described by Miiller, Klumpp and Lipowsky |3S].) If the trailing motor 
does not detach and it is assumed that the leading motor will reattach not too far from the position of 
the trailing head, then the system will proceed after having lost some ground. This sequence of events is 
common enough that Kunwar et al report that the stall force of a two motor system is subadditive when 
compared to the stall force of a single motor system. By contrast, our non-detachment model predicts a 
superadditive stall force. This is an important distinction to make because tug-of-war systems may spend 
substantial periods of time in a "stalemate" between forward and backward motors that are stalled, but are 
not exerting sufficient force to induce detachment |45j . 

5.3 Analysis including cargo dynamics 

Having developed some of the key two-motor results for the simplified setting neglecting cargo fiuctuations, 
we now proceed to extend them to include the cargo fiuctuations. To write a version of Theorem |5 . 1 1 that 
includes thermal fiuctuations of the cargo, we redefine the function G(r), which represents the instantaneous 
drift of a motor when the other is a signed distance r away. Rather than being written for a fixed cargo 
position, G now includes an average over the cargo position. We redefine the function G(r) as follows. From 
^ and ([T0|, 



G{r) EE 31 (x; 0) = 52 (x; 0) = A[g] ^) • (29) 
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Laser trap force (pN) 



Laser trap force (pN) 



Figure 5: Comparison of our stochastic averaging results with simulation of the effective force-velocity (left) 
and force-diffusivity (right) relationships of two cooperative motors attached to a common cargo. The direct 
stochastic simulations of the model equations ([T|) appear in blue with error bars denoting sample standard 
deviation. The stochastic averaging results, appearing as red solid curves, are computed from Eqs. (22 1 



and ( 23 1 for the velocity and diffusivity respectively, with G defined by ( 29 1 . The force- velocity curves 
are compared to the prediction of the force-balance theory {vg{6/2F^), magenta and dashed). The force 
diffusivity curves are compared to the diffusivity of the cargo when attached to one motor {(j'^/2, magenta 
and dashed). Parameter values are listed in Table 2.2 and we take the model ([3| for the instantaneous 
force- velocity relationship for a single motor. 



The calculation for the asymptotic average velocity V2{d) and diffusivity 'D2{0) proceeds exactly as in 
Section |5.1| with the effects of the cargo fluctuations appearing only in the redefinition of G. The accuracy 



of our stochastic averaging formula (Eqs. (22 1 and (29l) for the average velocity of two cooperative motors 
bound to a cargo is demonstrated in Figure [5| by its excellent agreement with direct numerical simulations 
of the underlying explicit stochastic dynamical model for the motor and cargo positions (Eq. [ij . We also 
contrast our theoretical force- velocity curve with what would be obtained by a simpler force-balance theory, 
in which the motors are always assumed to bear half the force applied to the cargo, taking into account neither 
cargo fluctuations nor the varying spatial configuration of the two motors | i45t i46j . We see that incorporating 
stochastic fluctuations of the cargo and relative motor positions does create a substantial change in the 
theoretical prediction of the effective transport of the motor-cargo system, even in the simplifled context of 
two identical cooperative motors without binding and unbinding dynamics (also noted for a discrete stepping 
model without cargo fluctuations in |75]) The effective diffusivity, seen in the panel on the right in Figure Isj 

2 — 

for the cargo attached to two motors is slightly more than half the value ^ of the effective diffusivity when 



attached to a single motor (111. The effective diffusivity increases somewhat as the applied force increases 
in either direction. 

Figure [6] compares the effective force- velocity relationship of the cargo when bound to one or two motors. 
In particular, we observe that the system with two motors is moving slower at low load forces than the 
single-motor-cargo system. On the other hand, at large applied forces, the two-motor-cargo system moves 
faster, and we see that for our parameter choices, the stall force for the two-motor-cargo system is about 
three times that of the single-motor-cargo system. These qualitative features have been previously noted in 
various models |37l 1721 [25] . usually through direct numerical simulations. Our analytical formulation allows 
us to provide more mathematical support to these observations. 

We begin by showing that G satisfies a concavity result similar to g. 

Lemma 5.3. Let g satisfy Assumption^ Then there exists an > such that for any s G (0, s^) we have 
that for all r > 

UG{r)+G{-r)) <G{0). (30) 
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Comparison of Velocity of 
One and Two Motor Systems 




Laser trap force (pN) 



Figure 6: Effective force-velocity curves for a single motor attached to a cargo (blue dashed, from Eq. ( |T3| ) 
and two cooperative motors attached to a common cargo (red solid, from (Eqs. (p2| and (29 1 ). Parameter 
values are listed in Table 2.2 and we take the model ([3| for the instantaneous force- velocity relationship for 
a single motor. 
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where we recall the function 77(2/) — {g{y) + g{~y)) from Assumption 1, which by hypothesis is decreasing 
for positive y (and therefore by even symmetry as a function of \y\). 

For each fixed and suitably small s, we produce a lower bound for the right-hand side in two cases: 1) 
rs < and 2) rs > f^,, where is the constant from Assumption [l] such that g"{y) < for all y < f^,. 

To address the first case, we rely on the concavity properties of g to deduce there exists an m > such 
that for y e [—/*,/*], ri"{y) < —to. Since y g [0, /*/2] implies for this first case that — /* < y — rs/2 < 
y + rs/2 < we obtain after a second order Taylor expansion: 



viy) 



rs 
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rs 
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forye [0,/,/2]. 



For all other y, we use the upper bound \ri{y) 
the sup norm over M. 



Uviy 



rs/2) + r,{y + rs/2))\ < ^|g"|oo, where 



Splitting the integral in (31 1 into two parts I = Ii + I2 with the integrals Ii and I2 being taken over 



[0, /*/2) and [/*/2, cx)) we have the estimate 



RHS (31) > 



8 



Zs > 



Zs > 



\9"\ 



where Zg is a normal random variable with mean and variance We can therefore choose an si > 

(independent of r) such that for all s < si, the right hand-side is positive. 

The Taylor approximation is not useful when r is not small, and so in the second case rs > we seek 
a uniform bound. Because 77 is a decreasing function of \y\, we see that for each fixed y, 77(7/ — a) + ri(y + a) 
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is a decreasing function of a for a > |?/|. Consequently, whenever < y < f* /2 and rs > f^,, we can deduce 
that r]{y) - {ri{y - rs/2) + ■q{y + rs/2))/2 > ■q{y) - {ri{y - /*) + ri{y + /*))/2 > if for some constant K >Q 
independent of r. 

Splitting the integral in ( pTj ) as in the first case, and combining the estimates from the previous paragraph 
with the global bound ri{y) — \{ri{y — rs/2) + 77(1/ + rs/2)) < 4|(7|oo, we obtain for the second case (rs > /*): 



RHS ([311 > > ^}W-4P{Z, > 



Once again, there exists an S2 > (independent of r) such that for all s < S2 the right-hand side is positive. 
Taking s.^, = min(si,S2) completes the proof. □ 

With this lemma we can now prove the following. 

Theorem 5.4. Let g satisfy Assumption^ Then there exist s^,,9^, > such that for any s < s* and \0\ < 6^,, 
we have V2(^) < Vi(^). 



Proof. We proceed as in the proof of Theorem 5.1 proving the property for 6 — and then extending it by 



continuity in 9 of an appropriate functional. Taking the computation of the average velocity from Section 



5.1 and applying (30) yields 



V2W = I I {G{r) + G{-r))7TR{r;Q)dr < I G{0)7TRir-0)dr = G(0) 



2 

Now, G(0) = .4[(7](0, ^), which we can rewrite as 



/ ,2, 1 /-oo 

^[5] (0, -j) = -^l [9i^y/^) + 9{-sy/V2) 



e ^ dy 



G(0) < ^ / (gisy) + g{-sy)) e~y' dy - 5(0) = Vi(0). 



□ 



Since g satisfies Assumption [T] (jm| , the integrand is decreasing and so 

poo 

V"" Jo 

5.4 Comparison to the high viscosity regime 

The results of the preceding sections establish that there is a clear interference effect when multiple motors are 
involved in cargo transport. However, the superadditivity of the stall force implies that when significant force 
is applied to the cargo, then multiple motors can cooperate and improve transport performance. It remains 
to discuss the performance of two motors in the presence of a high viscosity environment, in particular to 
address whether the forces encountered in vivo are sufficient so that two motors perform better than one. 
To this end, we show in Figure [T] the results of direct numerical simulations of the model (jlj for higher 
values of the cargo friction parameter 7 which might correspond to in vivo conditions. We see that the stall 
force for the two- motor-cargo system remains about the same in our model (> 20pN), about three times 
the stall force for a single motor. At low applied forces, the two-motor-cargo speed is even slower than that 
predicted by our asymptotic (small e) theory, and so is in particular slower than a single motor. Thus, the 
qualitative conclusions of our analytical theories for small cargo friction seem to remain true for the larger 
frictions expected in vivo. 

We do not yet have a direct analytical approach to identify the threshold viscosity that constitutes 
a boundary between the diffusion dominated and drag dominated regimes. But by comparing the direct 
numerical solutions to the the prediction of the stochastic averaging theory in Figure [7] we see that the 
unloaded cargo velocity begins to depart from the asymptotic theory when e ^ 10~^ (top panel), whereas 
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the overall structure of the force- velocity curve for a cargo attached to two motors shows significant deviations 
from the asymptotic theory only as e ^ 1 (lower panels). In either case, our stochastic averaging theory 
remains superior to the force balance theory. 

This is significant because when there is no external force applied and the drag force is 7 = 10"^ pN 
s/nm, we observe from Figures [t] that the velocities on one motor and two motor systems is roughly 275 and 
350 nm/s, respectively. Despite concerns that there may be interference between motors in vivo, it is more 
likely that the drag forces are sufficient that two motor transport is indeed faster. 




Figure 7: Top panel: Velocity of cargo attached to one or two motors as function of cargo friction parameter 
7. Results of direct simulation of ([ij are indicated by the blue solid curve (one motor) and magenta dashed 
curve (two motors), with error bars representing one sample standard deviation. The red dashed line (one 
motor) and black dot-dashed line (two motors) indicate the theoretical values corresponding to the e I 
limit, when the time scale of fluctuation of the cargo is much smaller than the time scale of the motor 
motion. The values of 7 presented correspond to a bead of size a = 500nm [18J and solvent viscosity ranging 
from the viscosity of water (at 300 K) to 1000 times that value (so e ranges from 3 x 10^"^ to 3). Lower 
panels: Numerical assessment of validity of stochastic averaging theory at higher solvent viscosities (left 
7 = 10~^pN • s/nm (77 = lO^ryo); right 7 — lO^^pN • s/nm (77 = 10'^77o) where r]o = 10~^pN • s/nm^ is 
the viscosity of pure water at room temperature). The meaning of the curves is the same as in Figure p] 



6 Conclusions 

We have developed a system of stochastic differential equations to describe the dynamics of multiple identical 
motors transporting intracellular cargo along a microtubule. We use a Langevin construction for the cargo 
equation and a diffusion approximation for the motor equations. 

We perform a dimensional analysis of our system of equations, and we find important dimensionless 
groups of parameters; this reveals that the dynamics in vitro are qualitatively distinct from what is likely 
to occur in vivo. In particular, the model suggests that in the experimental setting, the cargo is fluctuating 
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rapidly relative to the slower moving motors. 

In addition to the stochasticity found in the cargo dynamics, we study the impact of fluctuating motor 
positions. Using the randomness parameter calculated by we note that the motors' positions relative to 
each other are constantly changing. When multiple motors are engaged in transport of a single cargo, the 
force distribution among the motors is therefore state dependent. When the force- velocity relationship (Fig 
[2]) is nonlinear, it follows that the transport properties of teams of identical motors are non-trivial. We find 
that in the low viscosity environment, multiple motors are slower than one; however, for a fluid with viscosity 
similar to estimates of the cytoplasm |43] , drag force is sufficient that multiple motors perform better than 
a single motor. 

We also use our model to study the change in stall force as a function of motor copy number and viscosity 
of the fluid environment. This stall force is calculated to increase superadditively for multiple motors and 
also increases with the viscosity of the fluid environment. Neither of these results is in agreement with past 
work, and therefore a more careful study of motor dynamics near stall force is an important next step. The 
inclusion of binding and unbinding events will also clearly lead to a substantial change in the calculated stall 
force. In addition, we expect that to account for richer dynamics, the diffusion term in the motor SDE may 
need to be amended to a multiplicative, rather than additive, form. 
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A Extension to Nonlinear Spring Laws 

In the main text, we take a simple linear spring model ^ for the force law for the tails connecting the motors 
to the cargo. One could criticize this selection for a number of reasons. The tail actually has a complex 
hinged structure and would likely go "slack" under compression (as in the model of [31]). Also, we recall 
that in Eq. ([2]), r represents the separation between the motor and cargo projected along the microtubule, 
which is not the actual length of the tail in its real three-dimensional context. Finally, since the cargo is 
much larger (typically ~ 500 nm) than the motor or tether, the geometry of the motor's tail and location of 
its attachment to the cargo may be significant [31]. The effective one-dimensional force law may depend on 
how far away from the microtubule the connecting tail is bound to the cargo, and since this is potentially 
different for different motors, each motor's tail may have a different force law! The purpose of this appendix 
is to show that our modeling framework and stochastic averaging procedure can robustly incorporate these 
more realistic complications. 

The reason we have not pursued them in the main text, beyond the desire for a simple presentation 
of the key ideas, is that we could not identify a specific nonlinear force law that would be clearly bet- 
ter than the linear spring force law. The accuracy of experimental observations is too coarse to moti- 
vate a more detailed model for the connection; in fact controversy persists over what qualitative form 
the simple spring model should take [21]. Some authors prefer a wormlike chain model [25J F{r) = 
Kr + ((1 ~ kl/^c)^^ — 1) for |r| < that stiffens as the effective spring length approaches a critical 

contour length £ci beyond which the tether becomes taut and will not stretch any further. The framework 
we propose allows a general nonlinear spring structure to cover both of these, as well as other models for the 
force response of the tether attaching the motor and cargo. 
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To incorporate a more general tail force law, which could differ for each motor, the model equations 
can be modified ([T]) by simply replacing the tail force function F by a possibly motor-dependent tail force 
function Fi , and replacing the linear force law (Eq. |2| by a more general form: 



F{r) (r/Lc). 



(32) 



Here <i>i(^) describes an arbitrary nondimensionalized spring potential, Lc describes a length scale at which 
the spring force changes form, and F gives a characteristic magnitude of the tail force. The linear spring 
force law (|2| would correspond to a degenerate case in which $(^) — is not well defined (and 

cho 



could be chosen arbitrarily), and F — kL^. For a tail force law that is linear at small displacements 
but nonlinear after some length scale Lc (such as the wormlike chain model mentioned above), $ would 
be a more general function than quadratic, and F could be reasonably chosen as kLc, where k gives the 
constant of proportionality of the force against displacement in the linear regime. We will assume this type 
of model in the remaining discussion; other cases could be handled similarly with some modifications to the 
nondimensionalization procedure . 

Proceeding with nonlinear force laws that can be written in the form ( |32| (with F — kLc), and applying 
the same nondimensionalization as in Subsection [3] we obtain in place of Eq. ([t]) : 



N 



N 



di+dW,ii), 



(33) 



(34) 



with an additional nondimensional parameter defined as the ratio of the typical tail extension through 
thermal fluctuations to the length scale at which nonlinearity becomes important: 



A 



_ ^2kBT/K 



The dimensional analysis in Subsection [3] implicitly assumes this quantity is not large ( A < 1); otherwise 
using the length scale y/lk^T j k to characterize the typical tail length amounts to an inconsistent application 
of the linear spring force law beyond the regime of its validity. 

The stochastic averaging proceeds as before, provided that the $i(^) satisfy some technical growth and 
smoothness conditions (namely, $i should be and /p^exp{2 '^'^(v)dv\ dy — > ±oo as |a;| — >■ oo and 
exp |— 2 $^(w)d?;} dx < oo), to again yield the averaged equations Eq. (jsl, with the drift function for 
the motor coordinates and stationary distribution of the cargo modified as follows: 

g^(yi■,e)^ / 5 [sA"^$- (^(xi - z))] m^g(z)dz; 

N 



^xeW=^z'(X:^)exp<^-2 



^z + A-2^$, {X{x,~z)) 



exp < —2 



N 



^z + A-2^$, (A(x,-z)) 



The analysis in Sections |4.1| and |5.3| for developing the effective drift and diffusivity for a cargo bound 
to iV = 1 or iV = 2 motors proceeds in the same way, with just a few resulting changes in the formulas. 



The effective drift ( 13 ) of a single motor now involves a generally non-Gaussian average of the instantaneous 



force- velocity function: 

Vi{e)^Kz^ f 5[sA-i$'(A(a;-z))]exp|-2 9z + X'^<^> {X{x - z)) | dz, 



Kz = / exp <j -2 

/R 



(A(x-z))]} 



dz. 
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The formulas for the effective drift and diffusivity of a cargo bound to two motors are expressed in terms of 
the functions G±(r; 6), which would be redefined for nonlinear force laws as: 



G±{r;e) = Gi{r;9)±G2{r;e), 
gi{xi,X2;9) = Gi{xi - X2\0). 



References 



M. Asencio, G. Hooker, and H.O. Gao. Functional convolution models. 

R. Dean Astumian and Peter Hanggi. Brownian motors. Physics Today, pages 33-39, November 2002. 

G. Banting and S. J. Higgins, editors. Essays in Biochemistry: Molecular Motors, volume 35. Portland Press, 
London, 2000. 

Sebastien Bouzat and Fernando Falo. The influence of direct motor-motor interaction in models for cargo 
transport by a single team of motors. Physical Biology, 7:046009, 2010. 

O. Campas, Y. Kafri, K. B. Zeldovich, J. Casademunt, and J.-F. Joanny. Collective dynamics of interacting 

molecular motors. Phys. Rev. Lett, 97(3):038101, Jul 2006. 

N Carter and R Cross. Mechanics of the kinesin step. NATURE-LONDON-, Jan 2005. 

CM Coppin, JT Finer, JA Spudich, and RD Vale. Detection of sub-8-nm movements of kinesin by high- 
resolution optical-trap microscopy. Proceedings of the National Academy of Sciences of the United States of 
America, 93(5):1913, 1996. 

R.E. Lee DeVille and Eric Vanden-Eijnden. Regular gaits and optimal velocities for motor proteins. Biophysical 
Journal, 95:2681-2691, September 2008. 

Jonathan W. Driver, Arthur R. Rogers, D. Kenneth Jamison, Rahul K. Das, Anatoly B. Kolomcisky, and 
Michael R. Diehl. Coupling between motor proteins determines dynamic behaviors of motor protein assemblies. 
Phys. Chem. Chem. Phys., 12(35):10398-10405, 2010. 

Timothy C. Elston and Charles S. Peskin. The role of protein flexibility in molecular motor function: coupled 
diffusion in a tihed periodic potential. SIAM J. Appl. Math., 60(3):842-867 (electronic), 2000. 

ME Fisher and AB Kolomeisky. Simple mechanochemistry describes the dynamics of kinesin molecules. Pro- 
ceedings of the National Academy of Sciences, 98(14):7748, 2001. 

M. I. Freidlin and A. D. Wcntzcll. Random perturbations of dynamical systems, volume 260 of Grundlehren 
der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer- Verlag, New 
York, second edition, 1998. Translated from the 1979 Russian original by Joseph Sziics. 

John Fricks, Hongyun Wang, and Timothy C. Elston. A numerical algorithm for investigating the role of the 
motor-cargo linkage in molecular motor-driven transport. Journal of Theoretical Biology, 239(1) :33 - 48, 2006. 

J. Gagliano, M. Walb, B. Blaker, J.C. Macosko, and G. Holzwarth. Kinesin velocity increases with the number 
of motors pulling against viscoelastic drag. European biophysics journal, 39(5):801-813, 2010. 

Carla Goldman. A hopping mechanism for cargo transport by molecular motors on crowded microtubules. 
Journal of Statistical Physics, 140:1-15, 2010. 10.1007/sl0955-010-0037-2. 

Steven P. Gross, Michael Vershinin, and George T. Shubeita. Cargo transport: Two motors are sometimes better 
than one. Current Biology, 17(12):R478 - R486, 2007. 

Thomas Guerin, Jacques Prost, Pascal Martin, and Jean-Frangois Joanny. Coordination and collective properties 
of molecular motors: theory. Current Opinion in Cell Biology, 22(1):14 - 20, 2010. Cell structure and dynamics. 

N.R. Guydosh and S.M. Block. Backsteps induced by nucleotide analogs suggest the front head of kinesin is 

gated by strain, 2006. 

W.O. Hancock and J. Howard. Molecular Motors, chapter Kinesin: processivity and chemomechanical coupling, 
pages 243-269. Wiley- VCH, Weinheim, Germany 2003. 

P. Hanggi and F. Marchesoni. Artificial brownian motors: Controlling transport on the nanoscale. Rev. Mod. 
Phys., 81(l):387-442, January-March 2009. 



25 



[21] V. Hariharan and W.O. Hancock. Insights into the Mechanical Properties of the Kincsin Neck Linker Domain 
from Sequence Analysis and Molecular Dynamics Simulations. Cellular and Molecular Bioengmeering, 2(2):177- 
189, 2009. 

[22] AG Hendricks, E Perlson, JL Ross, HW Schroeder HI, M Tokito, and ELF Holzbaur. Motor coordination via a 
tug-of-war mechanism drives bidirectional vesicle transport. Current Biology, 2010. 

}23] N. Hirokawa. Kinesin and dynein superfamily proteins and the mechanism of organelle transport. Science, 
279(5350):519, 1998. 

}24] J. Howard, A. J. Hudspeth, and R. D. Vale. Movement of microtubules by single kinesin molecules. Nature, 
342 (6246): 154-8, 1989. 

125] Jonathon Howard. Mechanics of Motor Proteins and the Cyto skeleton. Sinauer Associates, Sunderland, MA, 
2001. 

126] J. Hughes, W.O. Hancock, and J. Pricks. Kinesins with extended neck linkers: A chemomechanical model for 
variable-length stepping. Bulletin of Mathematical Biology, 2011. 

127] J. Hughes, W.O. Hancock, and J. Pricks. A matrix computational approach to kinesin neck linker extension. 

Journal of Theoretical Biology, 269(1), 2011. 

128] D. Kenneth Jamison, Jonathan W. Drive, Arthur R. Rogers, Pamela E. Constantinou, and Michael R. Diehl. 
Two kinesins transport cargo primarily via the action of one motor: implications for intracellular transport. 

Biophysical Journal, 99:2967-2977, November 2010. 

129] R.Z. Khasminskii. On the principle of averaging in Ito's stochastic differential equations. Kybemetika (Prague), 
4:260-279, 1968. 

130] David Kinderlehrer and Michal Kowalczyk. Diffusion-mediated transport and the flashing ratchet. Arch. Ration. 

Mech. Anal, 161 (2): 149-179, 2002. 

[31] S Klumpp and R Lipowsky. Cooperative cargo transport by several molecular motors. Proceedings of the National 
Academy of Sciences of the United States of America, 102(48):17284, 2005. 

[32| Hiroaki Kojima, Etsuko Muto, Hideo Higuclii. and Toshio Yanagida. Mechanics of single kinesin molecules 
measured by optical trapping nanometry. Biophysical Journal, 73:2012-2022, oct 1997. 

[33] Anatoly B. Kolomeisky and Michael E. Fisher. Molecular motors: A theorist's perspective. Annual Review of 
Physical Chemistry, 58(l):675-695, 2007. 

134] CB Korn, S Klumpp, R Lipowsky, and US Schwarz. Stochastic simulations of cargo transport by processive 
molecular motors. The Journal of Chemical Physics, 131:245107, 2009. 

135] Peter R. Kramer, Juan C. Latorre, and Adnan A. Khan. Two coarse-graining studies of stochastic models in 
molecular biology. Communications in Mathematical Sciences, 8(2):481-517, 2010. to appear in Communications 
in Mathematical Sciences. 

[36] A Kunwar and A Mogilner. Robust transport by multiple motors with nonlinear force-velocity relations and 

stochastic load sharing. Physical Biology, 7:016012, 2010. 

137] A. Kunwar, M. Vershinin, J. Xu, and S.P. Gross. Stepping, strain gating, and an unexpected force-velocity curve 
for multiple-motor-based transport. Current Biology, 18(16):1173-1183, 2008. 

138] Y.A. Kutoyants. Statistical inference for ergodic diffusion processes. Springer Verlag, 2004. 

139] Adam G. Larson, Eric C. Landahl, and Sarah E. Rice. Mechanism of cooperative behaviour in systems of slow 
and fast molecular motors. Phys. Chem. Chem. Phys., 11:4890-4898, 2009. 

[40] C. C. Lin and L. A. Segel. Mathematics applied to deterministic problems in the natural sciences. Society for 
Industrial and Applied Mathematics (SIAM), Philadelphia, PA, second edition, 1988. With material on elasticity 
by G. H. Handelman, With a foreword by Robert E. O'Malley, Jr. 

141] Benjamin Lindner, Marcin Kostur, and Lutz Schimansky-Geier. Optimal diffusive transport in a tilted periodic 

potential. Fluctuation and Noise Letters, 1(1):R25-R39, 2001. 

142] Christian Maes and Maarten H. van Wieren. A Markov model for kinesin. J. Statist. Phys., 112(l-2):329-355, 
2003. 

143] CS Mitchell and RH Lee. A quantitative examination of the role of cargo-exerted forces in aoconal transport. 
Journal of theoretical biology, 257:430-437, 2009. 



26 



Alex Mogilncr, Timothy C. Elston, Hoiigyuii Wang, and George Oster. Molecular motors: Theory. In Christo- 
pher P. Fall, Eric S. Marland, John M. Wagner, and John J. Tyson, editors. Computational cell biology, volume 20 
of Interdisciplinary Applied Mathematics, chapter 12, pages xx+468. Springer- Ver lag. New York, 2002. 

M Miiller and S Klumpp. Tug-of-war as a cooperative mechanism for bidirectional cargo transport by molecular 
motors. Proceedings of the National Academy of Sciences, 105(12) :4609-4614, Jan 2008. 

M Miiller, S Klumpp, and R Lipowsky. Motility states of molecular motors engaged in a stochastic tug-of-war. 
Journal of Statistical Physics, Jan 2008. 

Melanic J.I. Miiller, Stefan Klumpp, and Reinhard Lipowsky. Bidirectional transport by molecular motors: 
Enhanced processivity and response to external forces. Biophysical Journal, 98(11):2610-2618, 2010. 

G Muthukrishnan, Y Zhang, S Shastry, and WO Hancock. The processivity of kinesin-2 motors suggests dimin- 
ished front-head gating. Current Biology, 19(5):442-447, 2009. 

Gayatri Muthukrishnan, Yangrong Zhang, Shankar Shastry, and William Hancock. The processivity of kinesin-2 
motors suggests diminished front-head gating. Current Biology, 19(5):442-447, Mar 2009. 

Igor V. Novozhilov. Fractional analysis. Birkhauser Boston Inc., Boston, MA, 1997. Methods of motion 
decomposition. Translation of the 1991 Russian original. 

G. A. Pavliotis. A multiscale approach to Brownian motors. Phys. Lett, A, 344:331-345, 2005. 

Thomas T. Perkins. Optical traps for a single molecule biophysics: a primer. Laser & Photonics Reviews, 
3(l-2):203-220, 2009. 

C Peskin and G Oster. Coordinated hydrolysis explains the mechanical behavior of kinesin. Biophysical journal, 
Jan 1995. 

C. S. Peskin, G. B. Ermentrout, and G. F. Oster. The correlation ratchet: A novel mechanism for generating 
directed motion by atp hydrolysis. In Van C. Mow, Farshid Guilak, Roger Tran-Son-Tay, and Robert M. 
Hochmuth, editors. Cell Mechanics and Cellular Engineering, pages 479-489. Springer- Verlag, New York, 1994. 

Filippo Posta, Maria R. D'Orsogna, and Tom Chou. Enhancement of cargo processivity by cooperating molecular 
motors. Phys. Chem. Chem. Phys., 11:4851-4860, 2009. 

Hong Qian. The mathematical theory of molecular motor movement and chemomechanical energy transduction. 
J. Math. Chem., 27(3):219-234, 2000. 

Peter Reimann. Brownian motors: noisy transport far from equilibrium. Physics Reports, 361:57, 2002. 

AR Rogers, JW Driver, PE Constantinou, DK Jamison, and MR Diehl. Negative interference dominates collective 
transport of kinesin motors in the absence of load. Physical Chemistry Chemical Physics, ll(24):4882-4889, 2009. 

M Schnitzer, K Visscher, and S Block. Force production by single kinesin motors. Nature cell biology, Jan 2000. 

S. Shastry and W. Hancock. Neck-Linker Length is a Critical Determinant of Kinesin Processivity. Biophysical 
Journal, 98:369, 2010. 

Shankar Shastrj' and William O Hancock. Neck linker length determines the degree of processivity in Kinesin-1 
and Kinesin-2 motors. Current Biology, 2010. 

Y Shtridelman, T Cahyuti, B Townsend, D DeWitt, and JC Macosko. Force-velocity curves of motor proteins 
cooperating in vivo. Cell Biochem Biophys, 52(l):19-29, 2008. 

Yuri Shtridelman, George M. Holzwaxth, Clayton T. Bauer, Natalie R. Gassman, David A. DeWitt, and Jed C. 
Macosko. In vivo multimotor force-velocity curves by tracking and sizing sub-diffraction limited vesicles. Cellular 
and Molecular Bioengineering, 2:190-199, jun 2009. 

G Shubeita, S Tran, J Xu, M Vershinin, S Cermelli, S Cotton, M Welte, and S Gross. Consequences of motor 
copy number on the intracellular transport of kinesin- 1-driven lipid droplets. Cell, 135:1098-1107, 2008. 

Anatoli V. Skorokhod, Prank C. Hoppensteadt, and Habib Salehi. Random perturbation methods with applications 

in science and engineering, volume 150 of Applied Mathematical Sciences. Springer- Verlag, New York, 2002. 

M Vershinin, BC Carter, DS Razafsky, SJ King, and SP Gross. Multiple-motor based transport and its regulation 
by Tau. Proceedings of the National Academy of Sciences, 104(1) :87, 2007. 

K Visscher, M Schnitzer, and S Block. Single kinesin molecules studied with a molecular force clamp. Nature, 
Jan 1999. 



27 



[68] Hongyun Wang. A new derivation of the randomness parameter. J. Math. Phys., 48(10):103301, 18, 2007. 

[69] Hongyun Wang. Several issues in modeling molecular motors. Journal of Computational and Theoretical 
Nanoscience, 5:2311-2345, 2008. 

[70] Hongyun Wang and Timothy C. Elston. Mathematical and computational methods for studying energy trans- 
duction in protein motors. J. Stat. Phys., 128(1-2) :35-76, 2007. 

[71] Hongyun Wang, Charles S. Peskin, and Timothy C. Elston. A robust numerical algorithm for studying biomolec- 
ular transport processes. J. Theoret. Biol., 221(4):491-511, 2003. 

[72] Ziqing Wang and Ming Li. Force-velocity relations for multiple-molecular-motor transport. Phys. Rev. E, 

80(4):041923, Oct 2009. 

[73] A. Yildiz, M. Tomishige, A. Gennerich, and R.D. Vale. Intramolecular strain coordinates kinesin stepping 
behavior along microtubules. Cell, 134(6):1030-1041, 2008. 

[74] Yunxiii Zhang and Michael E. Fisher. Dynamics of the tug-of-war model for cellular transport. Phys. Rev. E, 
82(1):011923, Jul 2010. 



28 



